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C^ ' Abstract 

In |[T], Sinopoli et al. analyze the problem of optimal estimation for linear Gaussian systems 

where packets containing observations are dropped according to an i.i.d. Bernoulli process, modeUng 

r \ a memoryless erasure channel. In this case the authors show that the Kalman Filter is still the optimal 

^^ I estimator, although boundedness of the error depends directly upon the channel arrival probability, p. 

f~| . In particular they also prove the existence of a critical value, pc, for such probability, below which the 

2 ' Kalman filter will diverge. The authors are not able to compute the actual value of this critical probability 

for general linear systems, but provide upper and lower bounds. They are able to show that for special 

cases, i.e. C invertible, such critical value coincides with the lower bound. This paper computes the 

J^ . value of the critical arrival probability, under minimally restrictive conditions on the matrices A and C. 
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A large wealth of applications demand wireless communication among small embedded de- 
vices. Wireless Sensor Network (WSN) technology provides the architectural paradigm to im- 



H I plement systems with a high degree of temporal and spatial granularity. Applications of sen- 

sor networks are becoming ubiquitous, ranging from environmental monitoring and control to 
building automation, surveillance and many others [i2J. Given their low power nature and the 
requirement of long lasting deployment, communication between devices is power constrained 
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and therefore limited in range and reliability. Changes in the environment, such as the simple 
relocation of a large metal object in a room or the presence of people, will inevitably affect the 
propagation properties of the wireless medium. Channels will be time-varying and unreliable. 
Spurred by this consideration, our effort concentrates on the design and analysis of estimation and 
control algorithms over unreliable networks. A substantial body of literature has been devoted 
to such issues in the past few years. In this paper we want to revisit the paper of Sinopoli et 
al. Q. In that paper, the authors analyze the problem of optimal state estimation for discrete- 
time linear Gaussian systems, under the assumption that observations are sent to the estimator 
via a memoryless erasure channel. This implies the existence of a non-unitary arrival probability 
associated with each packet. Consequently some observations will inevitably be lost. In this case 
although the Kalman Filter is still the optimal estimator, the boundedness of its error depends on 
the arrival probabilities of the observation packets. In particular the authors prove the existence 
of a critical arrival probability pc, below which the expectation of estimation error covariance 
matrix P^ of Kalman filter will diverge. The authors are not able to compute the actual value 
of this critical probability for general linear systems, but provide upper and lower bounds. They 
are able to show that for special cases such critical value coincides with the lower bound. 

A significant amount of research effort has been made toward finding the critical value. In [fTl], 
the author prove that the critical value coincides with the lower bound in a special case when 
the system observation matrix C is invertible. The condition is further weakened by Plarre and 
BuUo [3] to C only invertible on the observable subspace. In Q, the authors prove that if the 
eigenvalues of system A matrix have distinguished absolute values, then the lower bound is 
indeed the critical value. The authors also provide a counter example to show that in general 
the lower bound is not tight. 

Other variations of the original problem are also considered. In [[5]|, the authors introduce smart 
sensors, which send the local Kalman estimation instead of raw observation. In [[6]|, a similar 
scenario is discussed where the sensor sends a linear combination of the current and previous 
measurement. A Markovian packet dropping model is introduced in ^ and a stability criterion 
was given. In [8|, the authors study the case where the observation at each time splits into two 
parts, which are sent to the Kalman filter through two independent erasure channels. A much 
more general model, which considered packet drop, delay and quantization of measurements in 
the same time, is introduced by Xie and Shi H. 
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Another interesting direction to characterize the impact of lossy network on state estimation 
is to directly calculate the probability distribution of estimation error covariance matrix P^ 
instead of considering the boundedness of its expectation. In IfTOl . the author gives a closed-form 
expression for cumulative distribution function of P^ when the system satisfies non-overlapping 
conditions. In [[TTIl . the authors provide a numerical method to calculate the eigen-distribution 
of Pk under the assumption that the observation matrix C is random and time varying. 

In the meantime, lots of research effort has been made to design estimation and control 
schemes over lossy network, by leveraging the result obtained from above work. In 1121, the 
authors consider a stochastic sensor scheduling scheme, which randomly selected one sensor to 
transmit observation at each time. In [fT3l . the authors shows how to design the packet arrival 
rate to balance the state estimation error and energy cost of packet transmission. 

In a nutshell, we feel that derivation of critical value is not only important for analyzing 
the performance of the system in lossy networks, but also critical for network control protocol 
design. However, in a large proportion of the above work, the critical value is derived under the 
condition that C matrix is invertible or other similar conditions, which are not easy to satisfy 
for certain real application^. In this paper, we would like to characterize the critical value under 
more general conditions showing that it meets the lower bound in most cases. We also study 
some systems for which the lower bound is not tight and try to give some insights on why this 
is the case. 

The paper are organized in the following manner: Section |II] formulates the problem. Section Hill 
states all the important results of the paper, which will be proved later by Section HVl |Vl IVTl 
Finally Section IVIll concludes the paper. 

II. Problem Formulation 
Consider the following linear system 

Xk+i = Axk + Wk, 

(1) 
Uk = Cxk + Vk, 

where Xk G M" is the state vector, yk E W" is the output vector, Wk G M" and Vk E M™ are 

Gaussian random vectors with zero mean and covariance matrices Q > and R > 0, respectively. 

*C invertible implies that the number of sensors is no less than the number of states. 
May 17, 2010 DRAFT 



Assume that the initial state, Xq is also a Gaussian vector of mean xq and covariance matrix 
So > 0. Let Wi,Vi, Xf) to be mutually independent. Note that we assume the covariance matrices 
of Wi, Vi, xq to be strictly positive definite. Define |Ai| > IA2I > ■ ■ ■ > |A„| as the eigenvalues of 
A. 

Consider the case where observations are sent to the estimator via a memoryless erasure 
channel, where their arrival is modeled by a Bernoulli independent process {'jk}- According 
to this model, the measurement yk sent at time k reaches its destination if 7^,. = 1; it is lost 
otherwise. Let 7^ be independent of Wk,Vk,Xo, i.e. the communication channel is independent 
of both process and measurement noises and let P(7fc = 1) = p. 

The Kalman Filter equations for this system were derived in [[Q and take the following form: 

Xk\k = Xk\k-i + IkKkivk - Cxk\k-i), 

Pk\k = Pk\k~l — lkKkCPk\k~li 

where 

Xk+i\k = A.Xk\k, Pk+i\k = A.Pk\kA + Q, 
Kk = Pk\k-iC {CPk\k~iC + R) , 

2^0|-l = ^0) -Po|-i = ^0- 

In the hope to improve the legibility of the paper we will slightly abuse the notation, by 
substituting Pk\k-i with Pk. The equation for the error covariance of the one-step predictor 
is the following: 

Pfc+i = APkA^ + Q- ^kAPkC^iCPkC^ + RY^CPkA'. (2) 

If 7fcS are i.i.d. Bernoulli random variables, the following theorem holds ^: 
Theorem 1: If (A, Q^) is controllable, {C,A) is detectable, and A is unstable, then there 
exists a pc E [0,1) such thatElE 

sup EPk = +00 for < p < pc and 3Po > 0, (3) 

k 

EPk < Mp^ Vt for pe < p < 1 and VPq > 0, (4) 

^We use the notation supj, Ak — +00 when the sequence yl^ > is not bounded; i.e., there is no matrix M > such that 
Afc < M,Vfc. 

^Note that all the comparisons between matrices in this paper are in the sense of positive definite if without further notice 

May 17, 2010 DRAFT 



where Mp^^ > depends on the initial condition Pq > 0. 

For simplicity, we will say that EP^ is unbounded if sup;, EPf^ = +00 or EPf^ is bounded if 
there exists a uniform bound independent of k. 

III. Main Result 

In this section, we want to state all the important results for critical value, the proof of which 
can be found in later sections. Through out the rest of the paper, we always assume that the 
following conditions hold: 
(HI) {C,A) is detectable. 
{H2) A can be diagonalized. 
{H3) R, Q, So are strict positive definite. 

From Section [III it is clear that the critical value of a system should be a function of all 
system parameters, i.e. A, C, Q, R, Sq. However, the following theorem, the proof of which is 
in Section UVl states that the critical value does not depend on R, Q, Sq as long as they are all 
positive definite. 

Theorem 2: If i?, Q, Sq > are strictly positive definite, then the critical value of a system 
is just a function of A, C, and is independent of i?, Q, Eq. 

Since we have already assumed that R, Q, Sq are strictly positive definite, by Theorem |2l we 
can let R = Im, Q = In, So = In without loss of generality. Also since we assume that A can 
be diagonalized, we can always transform the system into its diagonal standard form. Hence, we 
assume that A is diagonal. We can also denote fc{A, C) as the critical value of system {A, C). 

When the dimension of A, C is large, which is often the case in reality, it is desirable to 
break the large system to several smaller blocks (or subsystems), which are easier to analyse. 
As a result, we define a block of the system in the following way: 

Definition 1: Consider the system {A, C) is in its diagonal standard form, which means A = 
diag{Xi, . . . , A„) and C = [Ci, . . . , C„]. A block of the system is defined as subsystem (Ax = 
diag{Xi^,...,Xi^),Cx = [C^,, . . . , C^J), I < ii < . . . < ii < n, where X = {ii,...,ii} C 
{1, . . . , n} is the index set. 

A special type of block, which we call equi-block, plays a central role in determining the critical 
value of the system and it is defined as 
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Definition 2: An equi-block is a block which satisfies |AjJ = . . . = |AiJ, and we denote it as 
{Ax^, CxJ, where Xg is the index set. 

Definition 3: V{A, C) is defined as the dimension of the largest equi-block of the system. 
The following theorem shows a basic inequality between the critical value of the original system 
and smaller blocks, which we will prove in Section HVl 

Theorem 3: Define fc{A, C) as the critical value for system {A, C). \f A = diag{\i, . . . , A„) 
diagonal and C = [Ci, . . . , C„], then 

fc{A,C)>MAx,Cx), (5) 

for all possible index set X C {1, . . . , n}. 

Before we continue on, we need to define the following terms: 

Definition 4: A system {A, C) is one step observable if C is full column rank. 

Definition 5: An equi-block is degenerate if it is not one step observable. It is non-degenerate 
otherwise. 

Definition 6: The system is non-degenerate if every equi-block of the system is non-degenerate. 
It is degenerate if there exists at least one degenerate equi-block. 

For example, if A = diag{2, —2) and C = [1,1], then the system is degenerate since it is an 

r 1 1 
equi-block and not one step observable. For A = diag{2, —2, 3, —3) and C = 

10 1 
two equi-blocks are {diag{2,—2), I2) and {diag{3,—3), I2) and both of them are one step 

observable. Thus, the system is non-degenerate. 

It can be seen that non-degeneracy is a stronger property than observability but much more 
weaker than one step observability. In fact, for a one step observable system, C matrix must have 
at least n rows, which implies yk is at least a M" vector. On the other hand, for non-degenerate 
system, the C matrix can only have X'(A, C) rows. In reality, V{A, C) is usually a small number 
comparing to n. 

In [[il, the authors proved that the critical value meets the lower bound when the system 
is one step observable. In this paper, we weaken the condition from one step observability to 
non-degeneracy. 

Theorem 4: If the system [H satisfies assumptions {B.\) — {H3) and the equiblocks of A 
associated to the unstable and critically stable eigenvalues are non-degenerate, the critical value 
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of the Kalman filter is 

Pc = max(l — |Ai|^^, 0) (6) 

where Ai is the dominant eigenvalue. 

For degenerate systems we can show that in general the critical value is larger than the one 
computed in theorem HI Nonetheless in this paper we will compute the critical value for second 
order degenerate systems. This includes a very practical case, involving complex conjugate 
eigenvalues. Let A = diag{\i, \2) G M?^"^. We can use the following theorem in conjunction 
with Theorem [3] as the building block to allow analysis of larger systems. 

Theorem 5: For a detectable system with A = diag{\i,\2), \\i\ > IA2I and R, Q, Sq > 0, 
the critical value is 

p, = MA,C) = max(l - |Ai|-2,0), (7) 

if the system is non-degenerate, or in other word, if one of the following conditions holds 

1) |Ai|>|A2|, 

2) rank{C) = 2. 

Otherwise the system is degenerate and its critical value is 

p, = /e(AC) =max(l - |Airi-°A/(-^-),0), (8) 

where Ai = A2exp(jV), and Dm{x) is the modified Dirichlet function defined as 

{0 for X irrational 
(9) 
l/q for X = r/q, r, g G Z and irreducible. 

IV. Properties of Critical Value 

In this section, we will prove Theorem [2] and [3l which demonstrate the relationship between 
critical value and system parameters. Throughout this section, we always assume that assumption 
(HI) - {H3) holds. 

First we want to prove the independence between critical value and the covariance matrix of 
the noise. 

Proof of Theorem \2\- Since R, Sq, Q > 0, we can find uniform upper and lower bounds 
a, a > 0, such that 

^Im < R < aim, ain < Sq < «/„, «/„ <Q< ain- 
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Let us define Pg = aJ„, Pq = ain, Pq = In, and 

Pk+i = AP^A^ + aIn - lkAP,C^{CP,C^ + aW'CP,A^, 

Pfc+i = AP,A^ + aIn - -fkAPkC^iCPkC^ + alJ-'CPkA^, 

P,%i = AP^A^ + h- iuAPlC^{CP;C^ + In.r^CPlA'. 

By induction, it is easy to check that P^ = aP^ and P^ = aP^ for all k. 
Also we know that Pq < Pq. By induction, suppose that P^, < P^, then 

£fc+i = ^Zfc^^ + «/n - 7fc^ZfcC'^(CPfcC^ + a/™)-'CP,A^ 

< APkA^ + «J„ - 7fcAPfeC^(CPfcC^ + «J„)"iCPfeA^ 

< APfcA^ + g - 7fcAPfcC^(CPfcC^ + P)-^CPfcA^ = Pfc+i. 
Hence, P;, < P^ for all k. By the same argument, P^ < P^ for all k, which implies that 

Pk = «n* < Pk < Pk = aP*. (10) 

Since a, a > 0, the boundedness of P^ is equivalent to the boundedness of P^. However by the 
definition of P^, we know that it is only a function of A, C which is independent of P, Q, Sq. 



We now want to prove that the critical value of a system is larger than the critical value of 
any of its blocks. 



>a 



Proof of Theorem \^ With out loss of generalitjO, we assume that X = {1, ...,/},/ < 
Let us define JT = {/ + 1, . . . , m} to be the complement index set of X. 
By Theorem [2l we suppose for the original system R = J^, Q = T,o = In- 
Let us define Pq = So = /„ and 

Pfc+i = APkA^ + 4 - ^kAPkC^iCPkC^ + l2m/2r'CPkA^, 
Cx 



771. 



where C 



C 



J 



p2mxn 



. Using Matrix Inversion Lemma, we can show that 



Pfe+l = A(P,-l+7fcC^C) 'A^ + Ir,. 



(11) 



fc+i 



A 



Pr' + 2^kC^C 



A^ + In. 



''if T — {1, . . . , m}, the proof is trivial. If I is an arbitrary subset of size I, we can always permute the states to make it 
equal {1, ■ ■ ■ ,1} 
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We know that Pq = Pq- By induction, suppose that Pk > Pk, then 

p.:' + ,.c^c - (P- . 2,,crc) ^ P- - P- .J\ '^ '^ 

\ ^2 '--'1 '-^2 ^"2 



CjCi 
CIC2 



= Pk' - Pk' + Ik 






= n"'-^fc"'-7JCi -G 



C? 



'T 



-C^ 



< 0. 



Hence, 



9 



Thus, by induction, Pk > Pk,'^k, which in turn proves 

fciAC)>fM,C). 



Now define Po,x = Ii, Po,j = Im-i and 

Pk+l,I 



AxPk,xAl + I1- ikAxPkxCl{CxPk.xCl + Im/2)-^CxPk,xAl. 



Pk+l,J — AjPk,jAj + Im-l — IkAj Pk,jC j{C jPk.jC J + Im/'2) CjPk^jAj. 



(12) 



(13) 



It is not hard to check that Pk 



Pk,x 
Pk,j 



, for all fc. As a result, Pk is bounded if and 



only if P/c X and Pk^j are both bounded. Combining (fT3l) . we know 

/c(A, C) > fM, C) = max{/,(Ax, Cx), fMj, Cj)}. 



V. Critical Value for Non-degenerate Systems 

This section is devoted to proving Theorem IH Before continuing, we would like to state 
several important intermediate results which are useful for proving the main theorem and have 
some theoretical value of their own. 

We will first deal with systems whose eigenvalues are all unstable. We will lift this restriction 
later in the paper. By "unstable" we mean that its absolute value is strictly greater than 1. We 
will call the eigenvalues on the unit circle critically stable and the ones with absolute values 
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strictly less than 1 stable. Since A is diagonalizable, we will restrict our analysis to systems with 
diagonal A. Also since some eigenvalues of A may be complex, we will use Hermitian instead 
of transpose in the remaining part of the article. 

Similarly to the observability Grammian, we find that the matrix Xli^i 'JiA~^^C^CA~^ plays 
an essential role in determining the boundedness of Kalman filter, which is characterized by the 
following theorem: 

Theorem 6: If a system satisfies assumptions (HI) — {H3) and all its eigenvalues are unstable, 
then supk EPk if finite if and only if E[(Xlj^i7i^"*^C'^Cyl~*)"^] exists, and it satisfies the 
following inequality 

oo oo 

aE[{J2 l^A-'^C^CA-')-^] < sup EPk < aE[{J2 l^A-'^C^CA-')-^], (14) 

where a,a>0 are constants. 

By manipulating Yl'i^i 'JiA~^^C^CA~\ we established the following result, which is essen- 
tially equivalent to Theorem IH but restricted on the systems whose eigenvalues are all unstable. 

Theorem 7: If a non-degenerate system satisfies assumptions {HI) — {H3) and all its eigen- 
values are unstable, then the critical value of the system is 

p, = i-iAiri 

If the arrival probability p > p^, then for all initial conditions, EP^ will be bounded for all k. 
Else if p < Pc, for some initial conditions, EPk is unbounded. 

Now we need to generalize this result to systems that have stable eigenvalues. The following 
theorem provides an important inequality for systems that have stable eigenvalues: 

Theorem 8: Consider a system satisfies assumption (HI) — {H3) with a diagonal A = 
diag{Ai,A2,A3), C = [Ci, C2, Csjfl. If Ai is the unstable part, A2 is the critically stable part 
and A3 is the stable part and CjS are of proper dimensions, then the critical value of the system 
satisfies the following inequality 

/e(A, C) < lim f,{dtagiaA,, a^a), [Ci, ^2])^. (15) 

a->-l+ 

Combining Theorem [3l |7] and [8l we will prove Theorem |4] in the last part of this section. 

'Note that there is no requirement of non-degeneracy for this theorem. 

*Note that diag{aAi,aA2) are unstable when a > 1. Hence the right hand side of the inequality can be computed by 
Theorem [7] 
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A. Proof of Theorem |^ 

In this subsection, we are going to prove Theorem [6l The key idea in the proof is to avoid 
analysing Riccati Equations, which were intensively studied in the previous works |[I1 ^, and 
try to formulate estimation error covariance Pk by maximum likelihood estimator. First let us 
write down the relation between 7j|/j and Xk'- 



IkVk 
7o2/o 

Xq 



IkCA-^ 
^-fc-i 



Xk+r 



IkVk 

Xo-a;o 



IkCA-^ 

loCA-^-^ 

^-k-i 











7oCA-2 ^^cA-' 



Wk 

Wi 
Wo 



(16) 



The rows where 7jS are zero can be deleted, since they do not provide any information to improve 
the estimation of Xk+i- To write (fT6l) in a more compact way, let us define the following quantities: 

A-^ ■■■ 



j^-k 
^-k-i 



A-^ 
A-2 A-^ 



nn{k+l)xn{k+l) 



Gk^e 



C 








c 



[n+m{k+l)]xn{k+l) 



Cfc — —GkFk 



Wk 

Wi 
Wo 



+ 



Vk 
Xo — Xq 



nn+'m{k+l) 



Tk 



CA~^ 

CA-^-^ 
^-fc-i 



g ]^(n+m(fc+l))xn y 



Vk 

Xo 



nn+m(k+l) 



(17) 



(18) 



(19) 



(20) 



May 17, 2010 



DRAFT 



12 

Define T^ as the matrix of all non zero rows of diag{'yklm, Ik-ilm, ■ ■ ■ , 7o-^m, In)- Thus T^ is 
a (m Yji=o li + '^) by (n + m(A; + 1)) matrix. Also define 

Yk = ^kYk, Tk = r^Tfe, Cfc = r^efc. 

Yfc, Tfc, Cfc are now stochastic matrices as they are functions of 7^, 7fc-i, . . . , 70. 
We can rewrite (fT6l) in a more compact form as 

Yfc = TfcXfc+i + efc. (21) 

From (I2TI) . we know that Y^ is Gaussian distributed with unknown mean TkXk+i and known 
covariance Cov(ek\^k)- Hence, we can prove the following lemmas: 

Lemma 1: If A is invertible, then the state estimation £^+11^ and estimation error covariance 
Pfc+i given by Kalman filter satisfy the following equations 

Xk+i\k = {f^Coviek\Tkr'fkr%''Cov{ek\Tkr'Yk, (22) 

Pk+i = {f^Cov{ek\TkVfk)-\ (23) 

Lemma 2: \f A = diag{\i, A2, ■ " " ) -^n). where |Ai| > IA2I > ■ ■ ■ > |A„| > 1, then F^F^ is 
bounded by 

P^^^/„(.+i) < FkF," < ^p^^-^4(.+i), (24) 

where F^ is defined in (flTI) . 

Lemma 3: If a system satisfies assumptions (HI) — {H3) and all its eigenvalues are unstable, 
then the error covariance matrix of Kalman Filter is bounded by 

«(Tf ffc)^^ < Pk+i < aif^nr' (25) 

where a,a gM are constants independent of 7^ and kU 

Proof of LemmaU} Given the observation Y^, by (|2TI) . we know that the maximum likelihood 
estimator of Xk+i is 

Xk+i\k = if^CovieklTky^fky^T^CovieklTky^Yk, 

^We abuse the notations of a, a, which will also be used several times later. These notations only means constant lower and 
upper bound, which are not necessarily the same in different theorems. 
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and the estimation error covariance is 



Pfc+i = {T^Cov{ek\Tkr'Tur'. 

Since Yk is Gaussian with unknown mean T^Xk+i and known covariance Cov{ek\Tk), we know 
that the maximum likelihood estimator is the optimal one in minimum mean error covariance 
sense. Thus, Xk+i\k and P^ given by maximum likelihood estimator are essentially the same as 
Xk+i\k and Pk of Kalman filter, which concludes the proof. ■ 

Proof of Lemma |2]- Notice that 

A 

^k 



-I A 



Therefore, 



H\~l 



A^ + I 


-A 




-A^ 








■ ■ . AA" + / 


-A 




~A^ 


AA^ 



(FkFl 



By Gershgorin's Circle Theorem [[T4l . we know that all the eigenvalues of (F^F^)"^ are located 
inside one of the following circles:|C — |Ajp — 1| = |Ai|, \C — |Ajp — 1| = 2|Ai|, |C — |Aip| = |Aj|, 
where ^s are the eigenvalues of (F^F^)^^. 

Since |Ai| > IA2I > ■ ■ ■ > |A„| > 1, for each eigenvalue of (F^F^)"^, the following holds: 



C > min{\Xi\^ + 1 - |Ai|, jA^p + 1 - 2|Ai|, |Ai| 



|A.|}, 



(26) 



and 



C < max{|Aip + 1 + |A,|, |A,p + 1 + 2|A,|, \X,\^ + \Xi\}. 



(27) 



Thus, < (|A„| - 1)2 < C < (|Ai| + 1)^ which in turn gives 



(fc+1)- 
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Proof of Lemma \3} Since Wi,Vj,Xo are mutually independent, 

Cov{ek) = Cov{GkFk[wk, ..., Wof) + Cov{[vk, • • • , t'l, Xq]) 

= GkFkdiagiQ, Q, . . . , Q)F^G^ + diag{R, i?, . . . , i?, Sq). 
Since we assume that Q = Im, -R = So = In, using Lemma [2l it is easy to show that 



;GkGf, + In+mk < Cov{ek) — GkFkF^ G^ + In+mk < 



— GfcG^, + In+mk- 



Since GfcGf = diag{CG^, ..., GG^ , GG"+Q, define n^ = A„,„(CC^) and ?!« = XmaxiGG") 
1 , we know that n^In+mk < G^Gj! < ncln+mk, which implies 

ar^r^ < Gov{ek\Tk) = TkGov{ek)Tl < aVkVl, 



where a 



(|Ai|+l)^ 



+ l,a 



>^G 



(|An|-l) 



jw + 1. Notice that 1^1;^ = /. Therefore, 



al < Gov{ek\Tk) < aL 
The above bound is independent of k and 7j, which proves 

aiT^TuY' < Pk+i = {T^Gov{ek\Tkr'ny^ < aif^'^ny'. 

Now we are ready to prove Theorem [6l 

Proof of Theorem 0- As a result of Lemma [3l we only need to show that 

1 



aE 



i=l 

Rewrite fjf f ^ as 



< sup-E 

k 



T^n 



<aE 



J2li^~'^C^CA~y^ 



i=l 



fc+1 



T^Tk = ^7fc+i_,A-^^C^CA-^ + A-('=+^)^A-(^^+i). 



(28) 



(29) 



i=l 



We know that Tl^Tu will have the same distribution as 



-k J-k 



fc+i 



Y,l^A^'''C''GA-' + A-('=+^)^A-(*^+i) 



i=l 



since 7jS are i.i.d. distributed. 

First we want to prove the right hand side of the inequality. Since all the eigenvalues of A 
are unstable, 

oo 

^A~'"G"GA~' <(3I, 



i=l 
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where /3 > is a constant. Thus 



i=l 



fc+1 



i=l V i=l 

fc+l oo 



i=l 

k+1 



i=k+2 

oo 



> Y^iA-'^C^CA-' + /3-^ Y^ -f,A-'^C^CA-' > mm{l,(]-^)J2l^A~'^C^CA- 



j=i 



=fc+2 



i=l 



Thus, 



E 



(T^f^Tfc)-^ <max(l,/3)E 



(^7,A-*^C^CA-T' 



i=l 



for all k. 



(30) 



which proves the right hand side of the inequality. 

Then we want proof the left hand side of the inequality. We know that 

k oo 



i=l 



i=l 



Thus, 



fc+i 



sup^[(Tf Tfc)-^] = supE[(^7iA-^^C^CA-' + A-(*^+i)^A-(fc+i))-i] 



i=l 

oo 



> sup E[{Yl^A-'''C''CA-' + yi-(fc+i)^A-{^^+i))-i] 



i=l 

oo 



= lim E[(y7,A-*^C^CA-* + A-('=+^)^A-('^+i))-i]. 

fc— >oo •^-— ' 

i=l 

Since A~''^A~'' = diag{\\i\~'^^, . . . , |A„|~^*^), A~''^A~'' is monotonically decreasing with re- 
spect to k in positive definite sense. Therefore (Xl^i Ti^'^^C^CA"* + A-('=+^)^A-('^'+i))-i is 
monotonically increasing. By the Monotone Convergence Theorem, we know that 

oo 

supE[(ffffc)-i] > lim E[{y" -f,A~'^C^CA-' + A-^''+^^^A~^''+^Y^] 

7, i-— VOO ' ^ 



1=1 

OO 



E[lim (V 7,A-*^C^CA-* + A-C^'+D^A-C^+i))-!] 



(31) 



1=1 



E[(5];7iA-*^c^cA-o-^] 



i=\ 



which proves the left hand side of the inequality. 
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B. Proof of Theorem 

In this part, we will manipulate (^^^ 7jA~*^C^CA~*)~^ to derive the critical value. The 
key idea is to use cofactors to find an upper bound of matrix inversion. For non-degenerate 
system, this upper bound will lead to a upper bound of critical value, which coincides with the 
lower bound in |[I1. 

Before we prove Theorem |71 we need the following lemmas. 

Lemma 4: For a non-degenerate system, it is possible to find a set of row vectors Li, L2, . . . , L„, 



such that LiC = [Ufi, 



,/i,„], where U^i 



1 and I 



«j 



if lA,- 



lA 



3\- 



Proof: It is simple to show that the lemma holds by using Gaussian Elimination for every 

equi-block. ■ 

Lemma 5: Consider that |Ai| > IA2I ■ ■ ■ > |A„|, U^i = 1 and lij = if i ^ j and |Aj| 



A.I- 



Define Azi = ii and Aij 



ij^iJ = 2,3, 



n. Then the determinant 



D 



l2,2^2 



'2,lA]^ /^ ^Ar, 



is asymptotic to YYk=i^k^^' ^■^■ 



ln,l^l " ln,2^2 



lim 



^l,nA„*^ 



I'n^n'^r. 



D 



. = 1. (32) 

Aii,Ai2,...,Ai„-^oo YYj^^^ A^**" 

Lemma 6: If a non-degenerate system satisfies (HI) — {H3) and all its eigenvalues are 
unstable, then following inequality holds 



lim sup 

Aii,...,Ai„— >-co 



Y:i=iA~'' C^CA~'J 



\2Aii 



</3In, 



(33) 



where /3 > is a constant, zi < Z2 < • • • < ^n G N, Azi = ii, Aij = ij — ij-i,j = 2, . . . ,n 

Proof of Lemma\5\ The determinant D has n\ terms, which have the form sgn{a) YYl^i 4,jfeA~' 
o" = (ji, J2, • • • , jn) is a permutation of set {l,2,...,n} and sgn{a) = ±1 is the sign of 
permutation. Rewrite (|32l) as 



*fc 



niLiA 






-«fc 



■jfc 



nLi K" 



^sgn{a)Wlk 



a k=l m=l ^ •'■■'■ 



n 



k=m ^jk 
\.k=m Afc 



k=l 

Aim 



(11^=1 Aj J 






(nLiA.)-^"---(nL.A.)-^^" 
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Now we can just analyze each term of the summation. Since |Ai| > ■ ■ ■ > |A„|, | nfc=m -^ifc I — 
I nfc=m'^fcl- P^^^*- consider that there exists some A;s such that |Ajj.| 7^ |Afc| and define k* to be 
the largest, which means lAj^., | ^ \Xk*\ and \\jj = |Afc| for all k greater than k*. Since \\k*\ is 
the smallest among |Ai|, . . . , |Afc|, we know that |Aj^, | > |Afc. |. Thus, 

1 lfc=fc* ^jk 



n 



k=k* ^k 



and 



lim 



n'>..n(ft 

k=l m=l ^ •'■•'■ 



k=m '^jk 



n 

k=m ^fe 



A, 



-Ai„ 



>1, 



n 

<in' 

k=l 



k,jk I 



lim 

Ai..*— >oo 



riLfc* '^. 



jk 



riLfc* ^k 



-Ai,. 



0. 



Then consider that if for all k, \\j^ \ = \Xk\, but (ji, . . . , j„) 7^ {1,2, ... ,n). Thus, there exists 
k* such that jk* 7^ A;*. Hence lk*,j^,., = 0. Therefore, these terms are always 0. 
The only term left is 

k=l m=\ VllA;=m ^ 

Thus, we can conclude that 



-Aj„ 



lim 

Aii,AJ2,---iAJn— ^00 



D 



nLi A, 



■«fe 



1. 



Proof of Lemma ^ Because the system is non-degenerate, by Lemma |4l we know that 
there exist Li, L2, ■ ■ ■ , L^, such that LjC = [/j,i, . . . , /i,„] is a row vector, /j j = 1 and /jj = if 
2 7^ j and I Aj I = \\j\. 
Define matrices 

^l,lAi ^ ^l,2A2 ^ ■ ■ ■ ^l,nA~*^ 



U 



^2,lAi ^ /2,2A2 ' 



^2,nA„*^ 



o = u-\ 



Define rT; = max(Amax(-^f -^1), • • • , \max{L^Ln)). Thus, Lf Lj < nilm, and 



i=i 



j=i 



nz 



n^ L 



j^~iiH^Hj^H . . . j^-inH qH j^H 



LiCA~'' 



L^CA-'- 



rii 



-U^'U, 



(34) 



(35) 
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and 



-1 



< ni {U"U) ^ = niOO" < nitrace{00^)ln 



^■'=1 ^ (36) 

= 71;^ Oij{0")j^iIn = ni'Y^ Oij X COnj{Oij)In = ^1^^ |Oijp/„, 

where conj{) means complex conjugation. 

Now by Lemma [51 we can compute the cofactor matrix of U and hence O = U~^. Define the 
minor Mij of U as the {n — I) x (n — I) matrix that results from deleting row i and column j. 
Thus 



0^, 



By Lemma |5j we know that 



lim 



det{U) 
det(f/) 

-ik 



1. 



Since Mij has the same structure with U, it is easy to show that 

n 

where pij is a constant. Thus, 



k=2 



lim sup 

AJi,...,Aj„— >oo 



RLiIAa 



|2AJ,, 



< lim sup 



niE.AO. 



i,j I ^,J 



AJl,...,Ain— >-CXD llfc=i l^fe 



2Aifc 



-/. 



lim sup rii V^ 



Ail ,...,A«„— >oo 



*J 



det(Mij^ 



det(f/) 



2 n 



fc=i 



.P^^^ 1 In 



<nAj2 



pIj 



«j 



nL2i'^fc 



«fe-ii 



n ^ -2fc 



nLiA 



/ni^^i'^"M"=^'E/'?/- 



fc=i 



«j 



By (|36]) . we can conclude that 



(3 = ni^p: 



2 
'id 



(37) 



(38) 



(39) 



Now we are ready to prove Theorem |71 

Proof of Theorem [71" By Lemma [6l we know that there exists ^i, . . . , ^n > 0, such that if 



A'ij > ^j, j = 1, . . . , n, then 



J2%iA-''^C^CA-'^ < 2(3Yl%i \>^j?'^"'- Define stopping 
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time ii = 'mi{i > ^1]% = 1}, ij = inf{i > ^j + ij-i\% = 1}, j = 2, . . . ,n. By the definition of 
ij and Lemma [6l it is simple to show 



-1 



J27^^""C^CA-'\ > Q]]A-*^^C^CA-^n >2(3Yl\Xj\^^'K (40) 

j=i J \j=i J j=i 

Therefore, E \{J2Zi%^~'^C^CA-'y^] is bounded if ^nLil^^'l is bounded. From 
the definition of random variable ij, we know that Aij are independent of each other. And 

P{Az, = k) = P(7.,.,+5, = 0, . . . ,7.,_i+5,+fc-i = 0,7.,_i+€,+fc = 1) = (l-p)'-^^-V, fc > e. 

(41) 
Now we can compute the expectation 






n 00 



(42) 



i=i fc=cj 

which is bounded if and only if 

|A,f(l-p)<l, j = l,2,...,n. 

We immediately know that the upper bound for the critical value is 1 — |Ai|~^. Combining with 
the lower bound given in [T], we can complete the proof. ■ 

Before we finish this subsection, we want to state the following corollaries about the estimation 
error covariance matrix and boundedness of higher moment of Pk- 

Corollary 1: If a non-degenerate system satisfies (HI) — {H3) and all its eigenvalues are 
unstable, then the estimation error of state Xk by using only observations yk-ij,yk-i2, • • • ? Z/fc-i„, 
where < ii <...<?„< /c and Aii = ii, Aij = ij — ij-i,j = 2, . . . , n, is bounded by 

n 

Cov{xk\yk-h,yk-i2, ■ ■ -^Vk-^iJ < /3 n IM'^^^'^n, (43) 

fc=i 

where /3 is a constant, provided that Aij are large enough. 

Corollary 2: If a non-degenerate system satisfies (HI) — {H3) and all its eigenvalues are 

unstable, then 

sup EP^ <oo^Pc>l- |Ai|~2^ 

k 

where g G N. 
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Remark 1: Note that in Corollary [H we do not assume any distribution of ii, ... ,in- Hence, 
this corollary allows us to take into account other packet drop models. 

C Proof of Theorem \S\ and |4] 

In this subsection, we will proof Theorem [8] and finally Theorem HI 

Proof of Theorem^ We consider the case where i?, Q, So are identity matrices. To prove 
the inequality, we will first show that fc{A, C) < fc{diag(Ai, A2), [Ci, C2]). Rewrite the system 
equations as 



Xk+1,1 
Xk+1,2 



Ai 



A, 



Xfc+1,3 = AsXk,: 



WhS, 



Vk 



Ci C2 



Xk,l 
Xk,2 

Xk,l 
Xk,2 



Wk,l 
Wk,2 



+ Vk + CsXk^s 



Now we want to build a linear filter. Since Kalman filter is the optimal linear filter, the critical 
value of Kalman filter should be no greater than our linear filter. 

Because ^3 is stable, we can just use £^,3 = A^Xk^i ^s an unbiased estimation of x^ 3 and the 
estimation error covariance is bounded. Now 0:^,3, Xfc_i^3, . . . ,Xo,3 become measurement noise 
and we know that 





{ 


Xk,3 


\ 




Gov 




Xk-l,3 




= 




K 


^0,3 


/ 





/ ■■ 


■ 4^"' 


A^ 


■■ 


I 


A2 


■■ 





I 





{ 


Wk-1,2 


\ 




" / ■■ 


^2 


Al 


Gov 




• 
















W^0,2 






■■ 


I 


A2 




\ 


^0,2 


/ 




. ■■ 





I 



(44) 
Using the same method in Lemma [2] and Theorem [3] , we can show that this covariance matrix 
is bounded by pl{k+i)n, where p is a constant independent of k. Thus, it is possible to find an 
i.i.d. sequence of Gaussian measurement noise v'f^ such that 



Gov 



(\ "' ' 


\ / 


<-l 


> Gov 


\l ^'0 . 


I \ 



G3Xk,3 
G3Xk-l,3 



\ I CsXq^s J / 



\ 
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We can build another system, 

Xk+1,1 

Xk+1,2 



Ai 



A, 



Vk 



Ci C2 



Xk,l 


+ 


Wk,l 


Xk,2 


J 


Wk,2 


Xk,l 


+ Vk + vl 


Xk,2 





(45) 



By the property of linear filter, the estimation error of Kalman filter for system (1451) will be 
greater than the one for the original system, which implies that 



fM,C) < fcidtag{A,,A2),[Ci,C2]). 



(46) 



Now define function 



g{X, A, 7) = AX A" + Q- -fAXC{CXC^ + Ry^CXA^, 

and a,/? to be scalars. Therefore 

g(X, aA, 7) - g{X, (3 A, 7) = (a^ _ /jS) [AXA^ - -fAXC^{CXC^ + R)-^CXA^]. 

Thus, g{X, a A, 7) is a non-decreasing function of a when a > 0. Since for system {A, C, R, Q, Sq), 
the error covariance matrix Pk follows recursive equation Pk+i = g{Pk,A,'jk)- By the mono- 
tonicity of g{X,aA,'jk), we know that Pk is also a non-decreasing function of a. Thus, the 
critical value for the system is also non-decreasing, which implies that 



fciAC)< lim fciaA,C). 



(47) 



The limit on the right hand side must exist because of the monotonicity of function fc. Combining 
(|46|) and (l47l) . we can finish the proof. ■ 

Proof of Theorem^- If the system does not have unstable and critically stable eigenvalues, 
then the proof is trivial. Otherwise by Theorem [3l we know that 



fciAC) > fMi,C,) = max{l - |Al|-^0}. 



By Theorem [8l 



fc{A,C) < lim fc{diag{aAi,aA2),[Ci C2]) = lim max{l — |q;Ai| , 0} = max{l — |Ai| ,0}. 

a— >1+ a-^l+ 

Hence, the critical value is 1 — |Ai|"^. ■ 
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VI. A Complete Characterization of Critical Value for Second Order 

Systems 

This section is devoted to a complete characterization of linear system with a diagonal 2 by 2 
A matrix. This kind of systems can be seen as the building blocks of larger systems and thanks 
to Theorem [3l we know that the critical value of such system will give a lower bound of critical 
value for larger system. To our surprise, the critical values of second order systems are in fact 
quite complex. 

For non-degenerate second order systems, we can directly apply Theorem H] to derive the 
critical value. Thus, we will focus on degenerate systems. Using the same strategy as the previous 
section, we will first deal with unstable degenerate system, then apply Theorem [8] to generalize 
the result to critically stable and stable systems. By the definition of degeneracy, we know that 
a detectable second order system is degenerate if and only if the following assumptions holds: 

1) A2 = Xiexp{j(p) , where j^ = —1 and (p E (0,27r). ((p ^ 0, otherwise the system is not 
detectable) 

2) rank{C) = 1. 

To simplify the notation, let us define A = |Ai| = IA2I, z = exp{j(p). The proof of critical 
value is divided into 2 parts. First we want to deal with the case when ip/2n is rational: 

Theorem 9: If a unstable second order degenerate system satisfies hypothesis (HI) — {H3) , 
then the critical value of the system is 



Pc 



l_|Ai| -29/9-1, (48) 



where (p/2iT = r/q , q > r and r,q eN are irreducible. 
Then we consider the case when ip/27i is irrational: 

Theorem 10: If a unstable second order degenerate system satisfies hypothesis (HI) — {H3), 
then the critical value of the system is 

p, = l-|Al|-^ (49) 

if Lp/27r is irrational. 

Proof of Theorem |3; By the properties of degeneracy, we know that rank{C^C) = 1. 
a^ ah 



Thus, C^C 



ah h'^ 



, where a, h are real constants. It can be also proved that a, 6 7^ due 
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to the detectability of {C,A). Since A = diag{Xi, A2) = diag{Xi, Xiz), we know that 



i=l 1=1 



a^\~^' ahX 


-2i^-i 








abX-^'z' b^X"^' 




/oo 


1 2" 


-i 




a 
b 



(50) 



Since a, 6 7^ 0, we know that E{Y.T=i liA-'^C^CA~')~^ is bounded if and only if 



E\Y.1.X- 



■2i 



i=l 



1 z-'' 

Z' 1 



< 00. 



Define 



E^-^-" 



j=i 



1 2-^ 
z' 1 



It is easy to show that 

00 
trace{Z) = 2^7iA-2\ 



1=1 

00 



det(S) = 5^7.A-^^ - E^^^"'^^^ X E^^^" 



-2^^^-^ 



A=l 



A=l 

00 00 



, j=l 

00 



00 00 



E^^^"" + 2E E 7.7.A--A-^-5:7A-^^-E E 7^7.^-"A- 



■ajV^i-i ^ ^i-n 



i=l 
00 00 



i=l j=i+l 



j=l 



4=1 J=j + 1 



5^ 5^ 7.7.A--A-^(2 

i=l j=i+l 



^*-J _ ^J-n 



(51) 



Define set Sg^oo = {^ G N|/ y^ kq,k E N} and S^^^ = {^ £ Sq,oo|/ < 0- Since z = exp{2rn/q) 
and q, r are irreducible, z-'* = 1 if and only if |j — z| ^ Sg,oo- It is easy to show that 

mi{2 - z'-^ - z^-'\\j - i\ €§,,00} = inf{2 - 2^ - 2-*|i = l,...,g-l} = 2-2cos(— ) > 0, 

and 

sup{2 - z'^^ - z^^'Wj -i\E §g,oo} = sup{2 -z' - z^'\i = 1, . . . , g - 1} < 4. 
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Thus, 

oo oo 

[2-2cos(27r/g)]^ J2 in,X~''X~'' < det{E) = J2 E 7.7, A"^^"^^' 

OO (^^^ 

<4E E 7.7.A"-A-^-. 

Define stopping time ri = inf{z G NI7.J = 1} and T2 = inf{j G N|j — n G Sg,oo,7j = !}• 

Thus, 

00 00 00 ^ 

A-^- < trace(S) = J] 7.A"^^ = E 7.A^^^ < ^ ^"'^ = I^^^"'" " ^^S) 

2=1 i=ri i=ri 

Now consider there exist two index a, b such that b > a,b — a e Sg,oo and 7a = 7b = 1- By 
the definition of ri, we know that ri < a. Suppose that b < T2, therefore ti < a < b < T2. By 
the definition of T2, a — Ti = kaQ, b — Ti = k^q. As a result, b — a= {h — ka)q, which contradicts 
with the fact b — a E Sg.oo- Therefore we can conclude that T2 < b. 

Thus, for all 7a7b = 1, 6 — a G Sg,oo,Ti < a,T2 < b, which gives 

00 00 

A-^-A-^-<X: E 7.7,A-^^A-^^- = E E 7.7,A-^^A-^^- 

i=l jf-ie§5,oo i=Ti jf>T2,j-iG§(?,oo 



:i-A- 

«=T1 J=T2 



Define ui, (T2 to be the eigenvalues of E. Thus, 

/^ In 1 1 (Ti + 0-2 traceiE) ^^^^ 

trace{E-') = a^' + a^' = ^ '- = 1/ . (55) 

(Ti(T2 det(^) 

By inequality (|52|) . (l53l) and (l54l) . it is easy to justify that Etrace{E~^) < cx) if and only if 

E ^ :^ J\ = EA^^^ = E(A2^iA2(-2-^i)) < 00. 

Now we need to compute the distribution of ti,T2 — ti. By definition, the event {ti = i} 
is equivalent to {71 = . . . = 7j_i = 0,7i = 1}. The event {t2 — ti = i], where i G Sg,oo, is 
equivalent to {7ri+j = 0, 7^-3 = 1}, for all j G §q,j. Since T2 — T1 only depends on 7ri+j, i G Sg^oo, 
^2 — Ti is independent of ri. The distributions of ti,T2 — ti are characterized by the following 
equations: 

P(ri = ^) = P(7i = . . . = 7,_i = 0, 7, = 1) = (1 - p)*- V, (56) 

and 

P(r2 - ri = z) = P(7,,+,- = 0,7.2 = 1) = (1 -p)"-lp, (57) 
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where j G Sg,i,i G §g,oo and |§qi| means the number of elements in S^^j. Thus, 

oo 

E (A^^^A'^"^-"^)) = E\^^' X ^A2(^2-^i) = ^j9(1 - pYX^' X J^ p(l - p)'^''''a2\ (58) 

The first series is a simple geometric series which is bounded if and only if p > 1 — 1/A^. 
Using the root test of convergence, X^ies P(^ — py^i^X^^ is bounded if and only if 



I / 2 * 

limsup '^9Wp(l -p)|s,,»lA2i = (1 -p)limsupA ^^ < 1. (59) 

|Sg,i|->-co |S5,i|->-oo 

Since |§q j| = \{i — l)(g — I)/*?], where \x~\ means the minimal integer that is no less that x, 
limsupig ._^oo| (i/|Sg,j|) = q/{q — 1). As a result, the second series convergences if and only if 

(1-p)A2^/(^^i) < 1, (60) 

which is equivalent to p > l — X~'^i/(i~'^\ Now we can conclude that the critical arrival probability 



IS 



p, = 1-A-^. (61) 



Proof of Theorem (TO]- The proof is quite similar to the proof of Theorem [9l The proof 
before (|5T1) still holds. However, in (|52|) we need to change the set §g,oo- Define set T^ oo = {I ^ 

N|2 - z^ - z~^ > e}, where e > 0. And T^^j = {/ G T^,oo|/ < «}• Therefore, dSll) becomes 

oo oo 

^E E 7^7,A-2^A-2^<det(H) = 5^^7.7,A-'*A-2^(2-^^-^-^^-0. (62) 

(1531 ), (1551 ) still hold if we change every set S^ ^ to T^ j. However only the left side inequality in 
(|54|) holds, because there is no guarantee that for all a,b satisfies b — a E T^oo^Ta = 76 = 1^ 
Ti < a,T2 < b always holds. Also in Inequality (|62|) . we only prove the left side inequality 
of (|52|) . As a result, EX'^'^'^ < oo will only be the sufficient condition for the boundedness of 
estimation error covariance. Following the rest of the proof, it can be derived that 

_2 i 

p > 1 — limsup A '^=''' (63) 

|T.,,|-+oo 

is sufficient for bounded estimation error. Since e can be any positive real number, we can 
conclude that 

Pc < 1 — lim lim sup A ^^^ . (64) 
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Now we only need to estimate i/|Tgi|. 2 — 2;* — z~^ = 2 — 2cos(i(/?). Thus 

2 — 2 cos{iip) > s -^ i(p ^ [2kiT — arccos(l — s/2), 2k7r + arccos(l — s/2)], A; G Z 
<^ i{ip/2'K) ^ [k- As,k + As], k eZ, 

where Ag = arccos(l — £:/2)/27r. 

Define A^^ = inf{i e N\i{ip/27r) e [k-2As, k+2As], k e Z}. Suppose that a(v?/27r), 6(^/27r), 6 > 
a, both belong to [A; - A^, A; + A^], A; G Z. Thus, (6 - a)(v9/27r) G [A; - 2A£, A; + 2AJ, A; G Z. 
By the definition of N^, we can conclude that b — a > N^, which implies that if a^ip /2tx) G 
[A; - Ag, A; + A^], then (a + l)(<^/27r), . . . , (a + iV^ - l)(v5/27r) ^ [A; - A^, A; + A^]. Therefore, 
if a ^ Tg^ooj then a + l,...,a + A^e — iG T^ oo- As a result, 

-— -^ > limsup -— - > 1. (65) 

Since Lp/2'K is irrational, lime_^o+ iV^ = 00. Therefore, 

1 
lim limsup— — - = 1. (66) 

By (l64l) and (l66l) . we can conclude that the critical arrival probability Pc satisfies 

Pc<i- A^^ 

which is exactly the lower bound in [[U. Therefore, we can conclude the proof. ■ 

Now we can proof the main theorem: 

Proof of Theorem \^ By Theorem |4l [9] and [TOl we know that the only case we need to 
prove is critically stable degenerate systems, which is trivial by directly applying Theorem [U ■ 

VIL Conclusions and Future Work 

In this paper we address the problem of state estimation for a discrete-time linear Gaussian 
system where observations are communicated to the estimator via a memoryless erasure channel. 
Following the work of Sinopoli et Al. [1], we were able to compute the value of the critical 
probability for a very general class of linear systems. The boundedness analysis in this paper 
can be easily generalized to general Markovian packet loss models and to the boundedness of 
higher moments of the error covariance. Future work will attempt at determining the complete 
statistics of the error covariance matrix of the Kalman Filter under Bernoulli losses. 
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